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ABSTRACT 



' Extended source size effects have been detected in photometric monitoring 

in ; 

t^J- of gravitational microlensing events. We study similar effects in the centroid 

o ' 
o 

OO 

' analytically. For a source with circularly symmetric limb-darkening profile, 



motion of an extended source lensed by a point mass. We show that the cen- 
troid motion of a source with uniform surface brightness can be obtained 



the centroid motion can be expressed as a one-dimensional integral, which 
can be evaluated numerically. We found that when the impact parameter is 
comparable to the source radius, the centroid motion is significantly mod- 



CZ3 

ified by the finite source size. In particular, when the impact parameter is 



smaller than the source radius, the trajectories become clover-leaf like. Such 
astrometric motions can be detected using space interferometers such as the 
Space Interferometry Mission. Such measurements offer exciting possibilities 
of determining stellar parameters such as stellar radius to excellent accuracy. 

Key words: astrometry - galaxy: halo - gravitational lensing - stars: fun- 
damental parameters 



1 INTRODUCTION 

Gravitational microlensing has produced many interesting scientific returns concerning the 
Galactic structure, dark matter and variable stars (see Paczyhski 1996a for a review). The 
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current (photometric) observations of microlensing events have one notable shortcoming: 
from the observed time scale of microlensing event one cannot infer the lens mass and dis- 
tance, since this measured quantity depends on the lens mass, distances to the lens and the 
source and the relative lens and source transverse motion. This degeneracy can be partially 
or even completely removed using some rare microlensing events (e.g., Gould 1996, 1997). 
A much more general way of lifting the degeneracy is conducting an accurate astrometric 
program together with the photometric observations ( Hosokawa et al. 1993, H0g, Novikov 
& Polnarev 1995, Miyamoto & Yoshi 1995, and Walker 1995); we refer to such an astromet- 
ric experiment as astrometric microlensing. Accurate astrometry can measure the angular 
Einstein radius and the relative parallax between the lens and source. With these two addi- 
tional constraints, the lens mass can be determined, even if we have no independent measures 
of the lens and source distances (Paczyhski 1998; Boden, Shao, & Van Buren 1998). This 
method also opens up the possibility of measuring the mass of nearby high proper-motion 
stars (Miralda-Escude 1996, Paczyhski 1996b). Boden et al. (1998) showed that the required 
astrometric accuracy can be reached by ground interferometers soon available on 8- 10m 
class telescopes such as KECK and VLT. The Space Interferometry mission (SIM), to be 
launched in 2004, offers the best hope for disentangling the lens degeneracy with its global 
astrometric accuracy of ~ /ias. Ground and space interferometers (for a review see Shao & 
Colavita 1992) therefore offers exciting possibilities of measuring the mass of lenses, which 
may turn out to be brown dwarfs, planets etc. 

Paczyhski (1998) pointed out that, using astrometric observations, it is possible to mea- 
sure the stellar radius to 1% accuracy for microlensed sources which show finite size effect. 
The photometric extended source effect has been studied by many workers (Witt & Mao 
1994; Gould 1994; Nemiroff & Wickramasinghe 1994). However, up to now, most studies of 
astrometric microlensing adopt a point source approximation; the only exception is Walker 
(1995), who briefly discussed the astrometric extended source effect. Since the photometric 
extended source effects have already been detected for several sources (95-BLG-30, Alcock 
et al. 1997; 97-BLG-28, Albrow et al. 1998 and 97-BLG-56) and in light of the possibility 
of measuring stellar radii and other stellar parameters to excellent accuracy, we study the 
astrometric extended source effects in detail. In §2, we first derive analytical expressions for 
the centroid motion for a source with uniform surface brightness. The centroid motion for 
sources with (realistic) limb-darkening profiles are then simplified into a one-dimensional 
integral, which can be computed easily numerically. In §3, we discuss briefly our results. 
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Throughout this paper, we ignore parallactic effects, since this subject has been covered 
in detail by other authors (e.g., Hosokawa et al. 1993, H0g et al. 1995, Miyamoto & Yoshi 
1995). 



2 RESULTS 

We shall work in angular coordinates and all angles are normalized to the Einstein ring 
angular scale given by (e.g., Schneider, Ehlers, & Falco 1992; Paczyhski 1996a) 
a na ( M \ 1/2 ( 10 kpc Y' 2 n _ D d D. 

6e = 0.9 mas ^_ j ' C =ZWV « 

where D S) D d , and M are the distance to the source and the distance to the lens (deflector), 
and the lens mass, respectively. For convenience, we first choose a coordinate system centered 
on the lens but later (for observation considerations) transform the results into a system 
where the source is at the origin (see eq. [15]). We use complex quantities to obtain a more 
compact representation of the lens equation (Witt 1990). In this formalism, the lens equation 
is simply 

C = *-| , (2) 

where z = x + iy are the coordinates in the deflector plane, ( = £ + irj the coordinates in the 
source plane, and z denotes the complex conjugate of z. Again, both z and ( are normalized 
to Oe given by Equation (1). (Note that for the case where z and ( are expressed in length 
units we obtain a different normalization in each plane (cf. Witt & Mao 1994)). 

Eq. (2) can be easily solved to yield two image positions, which are denoted by subscripts 
+ and — indicating their respective positive and negative parities (Witt & Mao 1994): 



c 

*+>- = 2 



l±Wl + 



ccj (3) 

The separation between the images are comparable to Oe when the distance between the 
lens and source is a few Einstein radii; for microlensing in the local group, D ~ 10 kpc, the 
separation is ~ mas, too small to be resolved. Therefore we can only observe the motion 
of the centroid of light for these two images. For a point source, the centroid motion forms 
an ellipse in a coordinate sytem where the unpreturbed source is at the origin (e.g., Walker 
1995; see eq. [16] below). In the next subsections, we first study the centroid motion for 
extended sources with uniform surface brightness and then for sources with more realistic 
profiles. 
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2.1 Analytical Results for a Uniform Source 

Let us now consider a circular source with radius r and constant surface brightness. The 
boundary of the circular source is simply described by = Co + re lip . Using eq. (3), we 
obtain a parametric representation for the image contours: 

Co + ret* 



1±W1 + 



9(<P) 



= x + -(<p)+iy+-((p), 0<v?<27r, (4) 

where we have assumed (without losing generality) that Co is on the positive real £-axis, and 
g(<p) = r 2 + Co + 2rCo cos <p. 

Since lensing conserves surface brightness, the magnification of a circular source with 
uniform surface brightness is just given by the ratio of the area of the images to the area of 
the source: 



2tt 



A*tot 



-y+m—, — + y-m 



dip 



dip 



dip 



(5) 



Note that the minus sign results from the fact that when ip changes from to 2n, the contour 
for the image with positive parity moves counter-clockwise whereas that with negative parity 
moves clockwise. This integral can be evaluated analytically; the results have been presented 
in Witt & Mao (1994). They are given below for completeness and comparison purposes. 
We can now calculate the centroid of light by weighting the image positions with mag- 



nification. For a circular source we obtain 
AO. 



— f 

nr 2 fi tot J 



dip 



(6) 



dip dip 
Obviously we have A0 y = because of the reflection symmetry with respect to the x-axis. 
After a partial integration the integral becomes 

/ (Co + rcosy?)(r + Cocosy?) 



AO 



3-n-rjj, 



9(<P) 



'1 



9(f) 



dip, 



(7) 



where g{<p) is defined under eq. (4). After some tedious algebra we find, similar to the case 
of magnification for a uniform source, the integral can be simplified into two expressions, 
one for the special case of Co = r and one for the more general case when ( ^ r. 
For Co 7^ r, one finds 
1 



AO, 



47rr 2 Co/Utot v /4 + (Co-r) 2 



7T 7T 7T 

ai_F( - , k) + a 2 E(-, k) + a 3 IT(-, n, k) 



(8) 
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(9) 



2nr 2 y/A + (Co - r) 2 <- 

where F, E, IT are respectively the first, second and third kind of elliptic integrals as defined 
in Gradshteyn & Ryzhik (1980), and 



2 

A0 X = (o = r, //tot = - 

7T 



n = . 4 ^° r , 9 and k = J ^ r= . (10) 

(Co + r) 2 V 4 + (Co - ?~) 2 

The coefficients are defined as follows 

0l = -(8 + r 2 + C 2 )(r + Co)(Co - r) 2 , ^ = -(8 - r 2 + C 2 )(Co - r), (11) 
a 2 = (4 + (r - Co) 2 )(r + Co)(r 2 + C 2 ), b 2 = (4 + (r - Co) 2 )(r + Co), (12) 
a 3 = 4(2r 2 C 2 + r 2 + C 2 )(r-Co) 2 /(r + Co), &3 = 4(1 + r 2 )(r - Co) 2 /(r + Co). (13) 
For Co = r , we obtain 

9 h 1 + r 2 

H — arctanr . (14) 

The centroid motion for this special case is particularly simple. Physically all images of the 
source boundary fall on two circles: one circle is just the Einstein ring centered at the origin 
and the other is centered on the source with radius (1 + Co) 1 ^ 2 - These two circles cross each 
other at two points, (0,1) and (0,-1). The contour for the positive parity images rotates 
counter-clockwise along the two right arc segments and that for the negative party images 
rotates clockwise along the two left arc segments. One can show the resulting central shift 
of these simple configurations is just Co- 
While it is possible to define the inertial system where the (unperturbed) source is at 
rest, e.g., by using global astrometry or the positions of nearby stars, it is often difficult 
or impossible to do so for the lens since it could be dark. It is therefore observationally 
convenient to express the centroid motion in a system where the (unperturbed) source is at 
the origin and the lens is moving, say, along the x-axis. Each trajectory is described by its 
impact parameter, p. The centroid position can be obtained as 

Ae' x = x-A6 x ^, AQ' y = v -AQ x ^- , (15) 
so (,o 

where the lens is at (x,p) and Co — {x 2 +p 2 ) 1 ^ 2 . Notice it follows from eq. (14) that for Co — r , 
the centroid shift in this coordinate system is zero. These results for extended sources should 
be compared with those for a point source (e.g., Walker 1995) 

= Ae » = -CfT2 ' (16) 
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When the lens moves from — oo to +00 parallel to the x-axis, the centroid motion forms an 
ellipse with the major axis along the x direction and the minor axis along the y-axis. 

The (topological) difference between the centroid motion of a point source and that of an 
extended source is approximately only a function of r/p. In Fig. 1, we illustrate the difference 
with one example for a source with r = 0.5. The centroid motions for four lens trajectories are 
shown with p = +3, +0.8, +0.5, and +0.1, respectively. When the impact parameter is much 
larger than the source radius (top left), the trajectory of an extended source nearly coincides 
with that for a point source, as expected intuitively. When the impact parameter decreases 
to ^ r (top right), the trajectory begins to deviate from the point source approximation, 
with the minor axis becoming smaller. For the special case when r = p (bottom left), the 
centroid motion goes through the origin; this occurs when the edge of the source touches 
the lens (cf. eqs. [14] and [15]). As the impact parameter is further reduced, the trajectory 
becomes clover-leaf like, going through the origin twice (bottom right) since the lens hits 
the edge of the source twice. 

2.2 Generalization to Sources with Center-to-Limb Variations 

The surface brightness of stars in reality shows center-to-limb variations. In this subsection, 
we extend our previous results (obtained for constant surface brightness) to such cases. To do 
this, we first obtain the magnification and centroid shift of an infinitely thin shell at radius 
R. We illustrate the derivation by considering the magnification of a shell. The image area 
of a source with radius R is given by A(R) = / u tot (i?) • rcR 2 , and that of a source with radius 
R + dR is given by A(R + dR) = /i to t(-R + dR) ■ ir(R + dR) 2 . The shell magnification is then 
just the differential image area divided by the shell area, i.e., [A(R + dR) — A(R)]/(2nRdR), 
which yields 



where the subscript s indicates the quantities are for a thin shell. Similarly for the centroid 
shift, we find 



Notice the extra magnification factors in A8 X arise because the centroid is weighted by light. 
For a source with circularly symmetric limb-darkening profile, I(R), the magnification and 
centroid shift are then simply 




(17) 




(18) 
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Figure 1. Four examples of astrometric trajectories showing finite source size effects. The source size is assumed to be 0.5 
Einstein radii (r = 0.5). For each example, the black dot marks the source position and the lens is moving from — oo to oo 
parallel to the x-axis. The centroid motion starts at the origin and moves counter-clockwise. The impact parameter (p) of the 
trajectory is shown at the left corner in each panel. For each example, the solid line shows the trajectory that takes into account 
the finite source size effect while the dashed lines shows that for a point source approximation. 

/ n s I(R)2nRdR J A6 Xia fj, s I(R)2irRdR 

/^limb = ^7 j A#2;,limb = . (19) 

jI{R)2nRdR J n s I(R)2irRdR 

o o 

These quantities can be easily evaluated for any limb-darkening profile and transformed into 
other coordinate systems, e.g., using eq. (15). 

Fig. 2 shows the predicted centroid motion for the first microlensing event (95-BLG- 
30) exhibiting (photometric) extended size effects discovered by the MACHO collaboration 
(Alcock et al. 1997). The source limb-darkening profile is modelled by 
I(R) I R* ( R 2 \ 

7 ^ = l- M i- W2 + Mi yi-^ + M2 ^l-^J, (20) 

with mi = 0.57, u 2 = 0.28 for the MACHO V band, and u x = 0.72, u 2 = 0.05 for the MACHO 
R band. The source radius and impact parameter are found to be r = 0.075, p = 0.055, 
whereas the Einstein radius angular scale is 9 E = 420/^as. For this event, the centroid 
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Figure 2. Simulated astrometric trajectories for the first microlensing event (95-BLG-30) that shows photometric extended 
source effects. The parameters used are taken from Alcock et al. (i997). The dashed ellipse is the centroid motion for a point 
source. The solid line shows the trajectory for a source with constant surface brightness, whereas the dotted and dot-dashed 
lines show the predictions for the MACHO V and R band, respectively. Note the scales are the two axes are different and 
should be compared with the astrometric accuracy of ~ fias of the Space Interferometry Mission. 



motion displays the clover-leaf like trajectory since we have p < r. The overall deviations 
can be easily detected with SIM with its ~ //as astrometric accuracy. On the other hand, 
the centroid motions for sources with different profiles are almost the same. We found that 
in almost all cases, the centroid motion for a source with limb-darkening profiles only show 
small differences from that for a uniform source. Such small differences will be difficult to 
detect astrometrically unless we have very dense time samplings. 



3 DISCUSSION 

Extended source effects in photometric monitoring of microlensing events have already been 
observed for several events. We studied the astrometric signatures of these events in future 
microlensing experiments using interferometry. For microlensing events that exhibit extended 
source effects photometrically, the impact parameter is expected to be comparable to the 
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source radius (e.g., Witt & Mao 1994); for such events, the centroid motion is significantly 
modified by the finite source size effect and the distinctive clover-leaf features can be easily 
detected with future interferometers such as SIM (cf. Figs. 1 and 2). These modifications 
should be taken into account when fitting the centroid motions^. This will add one more 
parameter to the fitting procedure, but this also means in principle the source size can 
be obtained independent of photometric observations. We conclude that using astrometric 
microlensing, it is possible to determine interesting stellar parameters such as stellar radius 
with ~ 1% accuracy (Paczyhski 1998). We are anxiously waiting for the arrival of astrometric 
microlensing. 
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